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1 Introduction 



Quantumchromodynamics (QCD) is the widely accepted theory to describe 
strong interactions between hadrons. This theory shows the well-known be- 
havior of asymptotic freedom. Furthermore, lattice calculation show a phase 
transition from the hadronic world to a system of free moving quarks and 
gluons, the Quark-Gluon-Plasma (QGP). Heavy-ion experiments at CERN- 
SPS, at recently started BNL-RHIC and at the yet to come CERN-LHC with 
energies of ^/s = 20, 200, 5500 A GeV respectively are designed to study the 
eventually formed QGP. But still there is lack of a dynamical description of 
both the transitions from hadrons to quarks and gluons and vice versa, derived 
from first principles from QCD. In this talk I present a classical, molecular- 
dynamical model, which contains explicitly the phenomenon of confinement 
and a dynamical mechajiiam for the formation of hadrons out of large system 
of quarks and gluons E'EfHu 



2 The chromodielectric model (CDM) 

start with the Lagrange-density originally imenled by Friedberg and Lee 
and intensively studied by several followers Ql3'BEJ. 

+ [d.aWa) U{a) . (1) 

The first term describes quarks, where g is a Dirac spinor with color-, spinor- 
and flavor-indices being suppressed and m is a mass-matrix for the different 
quark-flavor. = d'^ + igA^ is the covariant derivative, describing the 
minimal coupling to the gauge fields A^^ — k^A^^°- with coupling constant g 
and A", a = 1 ... 8, being the Gell-Mann matrices. 

The second term is the kinetic term for the gauge field in a medium, mediated 
via a dielectric function k((t). The color field tensor is given by = d^A'^ — 
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(9^^° where the /"'"^ are the structure constants of SU(3)c and 

one has -\F^^F^''"' = \ {E''-E'' - B'^-B"). E and B are the color-electric and 
-magnetic fields. 

The last two terms introduce a scalar field a with a quartic scalar potential 
U{a). This scalar field is designed to mimic the long-range behavior of non- 
perturbative QCD and is therefore purely classical. It acts like a mediuxn 
as in classical electrodynamics but with a dielectric constant k(o') < 1 M. 
The potential U{a) is adjusted to have a global minimum at the vacuum 
expectation value (VEV) a = ay and a local minimum at cr = 0. In the 
absence of color-fields, the scalar field takes on its VEV everywhere. 

2.1 Confinement mechanism 

The mechanism of confinement in the CDM is based on an interplay of the 
color-fields and the cr-field via the dielectric function k{(t) and the scalar 
potential which are shown schematically in fig. |l|. From the Lagrangian (|l|) 
one gets the equations of motion for the color-field. 

[D^,K{a)Fn=r , (2) 

where = g q^^ X"" qX"' / 2 is the color-current of the quarks.^ In an Abelian 
approximation the equations for the color-fields reduce to the usual Maxwell- 
equations dfj.{K{a)F^'^ = j'"'. The crucial point of the model is the choice 
of k{(t), which is supposed to contain all non-Abelian effects. It is unity in 
the absence of the cr-field and it vanishes when the scalar field takes on its 
VEV ay. If one considers a color-charge-distribution p with vanishing total 
color projections (in the Abelian directions 3 and 8), which we call a white 
cluster, a color-field is produced due to the Gauss-law W ■ {K{a)E"-) — p°- . 
The field is only allowed where k(ct) > (cr < cr„). To suppress the scalar 
field costs an energy [/(O) = B, and the vacuum exerts a pressure on the 
color-field. If the transition from K=ltoK — Oisa rapid one, then one is 
left with a well defined spatial region, where the scalar field nearly vanishes 
and the color field is non-zero. All color-field-lines start and end on charges 
inside this volume and therefore there are no Van-der-Waals-like interactions 
to other white clusters except for very short-ranged cr-efFects. For that reason, 
if eventually two white subclusters form inside the cluster, they can separate 
from each other. 

In addition, if the charge-distribution has a non-vanishing total charge, the 
field-energy deposited in this configuration is divergent, i. e. those configura- 
tion cannot be created. 

"Note that this current is not conserved, due to the color-carrying gluons. 
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2.2 Model equations 

As the confinement mechanism in our model depends solely on the specific 
choice of the dielectric function k((t) and the scalar potential U{a)^ we can 
neglect the spin-dependences in the quark-Lagrangian. Instead we replace it 
with a Lagrangian for classical, spinless charged particles that are coupled to 
the classical color-field. 



Lp = Y,'n^k\Jl-Sk Pn{x - Xk ) - (3) 



3^ 



" =9j2q^k<PN{x~x;k) (4) 



with the color-charge and the 4-velocity u'j^. In our numerical realization 
we deal with extended particles with a fixed Gaussian distribution 

p^(f-ffc) = (27rr^)-^/%-(-~---'=)V(2'-o^) , (5) 

where V < ar^ > = V^ro = 0.7 fm is chosen to fit the radius of the nuclcon. 
The color-current (^) is consistent with the originally derived equations of 
motion (|^) only in an Abelian sub-group of SU(3)c. This is related to the 
fact, that the gluons in QCD carry color as well, whereas tha-color-current 
carried by the gluons vanishes in the Abelian approximation ElO. 
To further simplify the numerical realization we neglect the magnetic fields 
_B". This is exact for static problems and for string like yoyo-excitations Ej. 
The two decoupled sets of Maxwell-equations reduce basically to the Gauss- 
law for each field. To summarize we now have the following equations of 
motion for the particles, the (electric) color-field and the (T-field 

Xk = (6) 



k = -qt J d\ (V0'^(f )) pn{x - Xk) (7) 
[K{a)^cj>'^{x;)) = -p'^ix) (8) 
U'{a)=^^a+\n'{a){yr{x))-{yr{x)) , (9) 



9V 



where the prime denotes differentiation with respect to a and <j)°'[x) is the 
electric potential which satisfies E"" = —'V(p°'{x), a G {3,8}. We choose for 
the scalar potential U{a) = B + aa'^ + ba^ + ca^ with B = (150MeV)'', a = 
(489.9MeV)2, b = -15901MeV, c = 163.1 and for the VEV = 61.1MeV. 
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The dielectric function is K{a) = ^cxp (^a{-^ — /3)^ + 1^ , where a = 7 and 
f3 = 0.4. The coupUng constant is chosen to reproduce the string-tension in 
an gg-configuration and takes on the value as = .g/(47r) = 2. 

2.3 Classification of particles 

In QCD the quarks and gluons are represented as triplet- and octet-states 
of SU(3)c respectively. In our classical simulation we assign classical charges 
g° to the quarks. These charges are the diagonal entries of the A", a S 3,8. 
Due to our approximation we neglect the non-Abelian part of the color-fields. 
Instead we treat these 6 gluon-fields as particles in the same formalism as the 
quarks except that these particle-gluons carry both a color and an anti-color. 
The corresponding charges for quarks, anti-quarks and gluons are depicted in 
fig. 1^. E.g. the particle r in that scheme is a quark of color red and color- 
charges 1 and 1 / -\/3 with respect to the 3- and the 8-field respectively. 



As we have mentioned in section 2.1, the dynamics of our model forces the 
charged particles into white clusters and the separation into white subclusters 
is allowed. We regard as hadrons only white clusters which cannot be divided 
into smaller ones. It turns out, that there is only a finite set of those irreducible 
white Clusters (IWC). The IWCs consist either of a quark and an anti-quark or 
of three quarks (anti-quarks), both with some gluon admixture, or they consist 
only of gluons and can therefore easily be interpreted as mesons, baryons and 
glue-balls respectively. 

For the particle masses we use constituent quark masses to fit the low-lying 
hadronic spectrum. As our model does not depend on isospin, we treat the u- 
and the d-quark as degenerate particles with the same mass. Only the light 
pion, which is assumed to be a Goldstone-boson of chiral symmetry breaking, 
does not fit in our constituent-mass-scheme and thus we do not incorporate 
pions in our model. The quark masses are fixed to be mu,d = 400MeV, 
nig — 550MeV and rric — ISOOMeV. Thus the masses of the lowest hadronic 
states are simply the sum of the quark masses of the IWC. For the gluons 
we take a mass nig = 700MeV to reproduce the lightest glueball-mass of 
1400MeV. 

3 Hadronization 

To simulate the hadronization out of a QGP, we start with an ensemble 
of quarks and massive gluons, distributed homogeneously in a sphere of 
radius R = 4fm in real space and according to a Boltzmann-distribution 
with initial temperature Tq = 160MeV in momemtum-space. The rela- 
tive number of different particles is given through the distribution Ni oc 
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di exp(— TOi/To) ( "2^° )^^^' where di is a (spin, isospin and color) degener- 
ation factor. The colors are chosen randomly with the constraint of overall 
color-neutrality. After solving the Gauss-law for the color-fields in the first 
time step the system is driven by the equations of motion Due to the 
initial momenta the particles tend to leave the system but are bound due to 
the formation of color-strings. In this way the particles are reorganized to 
form in a first step white clusters and finally only IWCs. This scenario is 
shown in fig. |3[ The hadrons have invariant masses M?^^ — + with 
particle energies Ei and momenta pi 



E 



E, 



+ / d^x (ia^ + i(V(7)2 + U{a) + \K{a)E"'E''^ (10) 



P = 



E 



Pi 



(11) 



The resulting mass distribution is shown in fig. The curves are fits to 
a Hagedorn distribution dn/dm (x m""'"'^/^ exp(— TO(l/r — l/T/i)) and to an 
inverse power-law-distribution dn/dm oc . If we assume Th — 160MeV 
and a = —3 or a = —3/2 in the Hagedorn case, we get a hadronic temperature 
T = 146MeV and T = 125MeV respectively. 





Figure 1. The scalar poten- 
tial U (cr) and the dielectric 
function k((t). 
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Figure 2. The color scheme in the Abelian approximation. 
On the 3(8)-axis is given the particle charge with respect 
to the 3(8)-field. The particle on the two triangles (r, g, 
b) and (f , g, b) represents the quarks and anti-quarks and 
those on the hexagon the 6 massive gluons. 
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Figure 3. Different timesteps of the hadronization process. The simulation starts with a 
droplet of quarks and gluons and ends after about 200 fm left with IWCs only. The particle 
labels desribe the flavor content and grey shaded particles are already detected IWCs. 
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Figure 4. The mass distribution of the resulting hadrons. The curves are fits to a 
Hagedorn distribution with a = —3/2 and a = —3 and b = (1/T — l/T^) and to 
an inverse power-law-distribution. 
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